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Abstract 

The  angular  distribution  of  radiation  from  a  rod-pinch 
diode  is  examined  with  coupled  PIC-Monte-Carlo 
simulations  using  the  LSP  code[l]  and  compared  with 
angular  dose  measurements  on  the  ASTERIX  pulsed- 
power  generator.  For  the  4-  to  6-MV  voltage  range 
accessible  in  the  experiments,  both  the  measured  and 
simulated  doses  at  1  m  were  25  to  38  rad-LiF  in  the 
forward  direction  (0°).  The  calculated  dose  at  90°  is  1.2- 
1.8  times  larger  than  the  forward-directed  dose  depending 
on  both  the  voltage  and  rod  diameter.  This  is  also  in 
agreement  with  the  measured  angular  variation[2].  At 
6  MV,  the  simulations  show  that  the  dose  near  1 80°  is 
more  than  twice  the  forward-directed  dose.  At  these  high 
voltages,  electrons  approach  the  anode  rod  primarily  at 
angles  close  to  180°  resulting  in  peak  radiation  in  the 
backward  direction.  This  suggests  that  negative  polarity 
should  be  considered  for  voltages  at  or  above  6  MV  to 
maximize  the  extracted  dose.  Simulations  show  that  the 
0°  dose  is  higher  with  tapered  rods  and  increases  with  rod 
diameter.  Simulations  at  1 0  MV  show  that  the  dose  near 
180°  is  more  than  ten  times  higher  than  the  0°  dose.  At 
10  MV,  the  dose  near  180°  at  1  m  from  a  2-mm  diameter, 
blunt  rod  is  calculated  to  be  -800  Rad  (FiF). 


I.  RADIATION  MEASUREMENTS  AND 
CALCULATIONS 

A  standard  rod-pinch  diode  consists  of  a  blunt  or 
tapered  tungsten  or  gold  anode  rod  that  protrudes  a  few 
cm  beyond  an  annular  cathode  as  shown  in  Fig.  1.  For  the 
Asterix  experiments[2,3]  analyzed  here,  the  anode  was 


either  a  1-  or  2-mm  diameter  tungsten  or  gold  rod  that 
extended  1.6  cm  beyond  the  cathode  annulus  and  was 
tapered  to  a  point  over  the  last  1.0  cm.  The  electron 
positions  from  a  particle-in-cell  (PIC)  simulation  (see 
Fig.  1)  show  that,  because  of  the  strong  self-magnetic 
field  and  the  emission  of  protons  from  the  anode  (not 
shown),  electrons  are  strongly  focused  onto  the  tip  of  a 
high-atomic-number  anode  rod.  The  result  is  an  intense 
source  of  x-rays  with  a  source  diameter  comparable  to  the 
diameter  of  the  rod  when  viewed  along  the  rod  axis.  [4] 


Figure  1.  Geometry  of  the  rod-pinch  diode  with  electron 
positions  from  a  PIC  simulation. 

The  angular  distribution  of  the  dose  is  measured  at 
several  angles  using  both  thermoluminescent  dosimeters 
(TFDs)  and  Si  pin  diodes.  Time  histories  of  the  radiation 
are  measured  with  an  array  of  six  lead-collimated  Si  pin 
diode  detectors.  The  angular  locations  of  the  detectors  are 
shown  in  Fig.  2.  These  radiation  diagnostics  are 
described  in  detail  in  Ref.  2.  The  tip  of  the  rod  is 
positioned  near  the  center  of  a  10-cm-radius,  1-cm-thick, 
hemispherical  aluminum  shell  so  that  all  radiation 
measurements  are  attenuated  through  1  cm  of  aluminum. 
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Figure  2.  A  schematic  illustrating  the  positions  of  the 
radiation  diagnostics  used  in  the  Asterix  experiments. 


The  measurements  show  that  the  dose  increases  with 
angle  and  the  pin  diodes  at  large  angles  are  subject  to 
saturation  effects. [2]  Therefore,  the  pulse-shape  from  the 
pin  diode  located  at  -20°  is  used  to  compare  the  measured 
and  calculated  dose-rates.  There  is  no  evidence  of 
saturation  effects  in  the  smaller  signals  from  this  detector. 
The  time-histories  of  the  radiation  at  the  various  TLD 
locations  is  then  inferred  by  scaling  the  -20°  pin  diode 
pulse-shape  so  that  the  time-integral  is  the  TLD  dose  at 
that  location.  This  scaling  is  appropriate  because  the 
measured  TLD  doses  were  observed  to  scale  linearly  with 
the  dose  from  the  pin  diodes  in  a  previous  experiment  on 
Asterix. 

The  angular  distribution  of  the  radiation  from  a  rod- 
pinch  diode  is  calculated  with  the  particle-in-cell  (PIC) 
code,  LSP[1,5].  LSP  incorporates  electron  scattering, 
energy  loss,  and  bremsstrahlung  production  algorithms 
from  the  Integrated  Tiger  Series [6]  (ITS)  of  electron- 
photon  Monte-Carlo  transport  codes.  The  method  which 
is  used  to  calculate  the  angular  dose  distribution  is  similar 
to  that  used  to  calculate  the  forward-directed  (0°)  dose  for 
voltages  up  to  4  MV. [5]  First,  LSP  is  run  to  steady  state 
for  a  number  of  voltages.  Once  a  steady-state  voltage  is 
achieved,  a  list  of  electron  positions  and  momenta 
incident  on  the  surface  of  the  rod  is  generated.  The 
electron  list  is  post-processed  with  the  ITS  code, 
CYLTRAN,  to  calculate  the  angular  distribution  of  the 
radiation  spectrum  emerging  from  the  rod.  An  alternative 
method  of  calculating  the  radiation  spectrum  utilizes  the 
capability  of  LSP  to  generate  a  list  of  photons  that 
contains  their  energy,  direction,  and  location  in  the  rod  as 
they  are  created.  The  angular  distribution  of  the  radiation 
spectrum  is  then  calculated  by  using  CYLTRAN  to  model 
self-absorption  in  the  rod.  This  is  the  method  used  in 
Ref.  5.  The  two  methods  provide  an  excellent  cross¬ 
check  and  we  have  found  that  both  methods  produce 
essentially  the  same  results. 

The  angular-dose  distribution  at  1  meter  is  calculated 
from  the  radiation  spectrum  using 

Dosem  =  r R(hv)F(hv)^.(W  d(M  ,  (i) 

Qe  Jo  d(hv) 


where  Qe  is  the  incident  electron  charge  on  the  rod,  R(h  v) 
is  the  spectral  response  of  the  detector  [dose  per  photon 
flux],  F(hv)  is  the  attenuation  from  any  filters  between  the 
source  and  the  detector,  and  dOs/d(hv)  is  the  angular 
distribution  of  the  photon  flux  spectrum  per  incident 
electron  charge  calculated  from  the  LSP/CYLTRAN 
model.  Both  the  detector  response  and  the  filter 
attenuation  are  calculated  using  tabulated 
attenuation/absorption  coefficients. [7]  This  approach 
avoids  the  problem  of  poor  photon  counting  statistics 
associated  with  the  small  solid  angle  subtended  by  the 
detectors  incurred  with  direct  modeling  of  the  detectors 
with  ITS.  However,  Eq.  (1)  omits  radiation  scattered  in 
the  filters  that  can  subsequently  reach  the  detector  and  it 
can  overestimate  the  energy  absorbed  energy  in  the 
detector  if  it  is  not  properly  equilibrated.  Both  of  these 
effects  are  small  in  this  experiment  and  the  forward  dose 
predicted  by  Eq.  (1)  agrees  with  that  calculated  in  Ref.  5 
where  the  detectors  are  modeled  directly  with  ITS. 

The  incident  electron  charge  on  the  rod  in  Eq.  (1)  must 
be  calculated  in  terms  of  measurable  quantities  to  be 
useful  in  comparing  with  experiments.  This  charge  is 
determined  from  the  measurables  as 

Qe  ~  Qt  N passes  (1~^ior/^r)  ?  (2) 

where  Qt  is  the  total  absorbed  charge  (i.e.  the  integral  of 
the  current)  in  the  diode  (both  electrons  and  ions),  Npasses 
is  the  number  of  passes  the  electrons  make  in  the  rod 
defined  by  the  ratio  of  incident  electron  charge  to 
absorbed  electron  charge,  and  Ij0n/lT  is  the  ion-current 
fraction.  Npasses  can  be  obtained  from  either  CYLTRAN 
or  LSP  and,  for  1-  and  2-mm  diameter  tapered  rods,  the 
simulations  show  that  electrons  make  between  2  and  10 
passes  through  the  rod  as  the  voltage  increases  from  2  to 
10  MV.  The  ion  current  fraction  from  LSP  varies  from 
25%  to  45%  when  protons  are  the  ion  species  and  it  is  a 
weak  function  of  the  voltage.  For  blunt  rods,  electrons 
make  only  2  to  3  passes  through  the  rod  and  the  ion 
current  fraction  is  reduced  to  25%  to  35%  for  voltages 
between  2  and  10  MV  since  electrons  spend  less  time  in 
the  diode.  To  apply  Eq.  (1)  to  time-dependent 
experiments,  the  dose-per-charge  is  interpreted  as  the 
dose-rate  divided  by  the  total  current  so  that 

Dose/QT=(Dose-Rate)/IT ,  (3) 

where  Ij  is  the  total  current  in  the  diode.  This  relationship 
is  valid  since  the  LSP  calculations  are  at  steady-state. 

Results  from  angular  dose  calculations  for  a  2-mm 
diameter  tapered  rod  are  shown  in  Fig.  3.  The 
calculations  show  that  the  angular  variation  of  the  dose- 
rate  efficiency  ( Dose-Rate/Ij )  is  small  for  voltages 
below  4  MV  but  increases  rapidly  above  4  MV.  The 
simulations  also  show  that  the  0°  dose-rate  efficiency 
becomes  insensitive  to  changes  in  the  voltage  for  voltages 
above  6  MV.  This  insensitivity  is  due  to  electron  angles 
of  incidence  on  the  rod  approaching  180°  as  the  voltage 
increases  which  increases  the  x-ray  emission  at  large 
angles  increases  relative  to  the  smaller  angles.  The 
model  predicts  that  the  dose-rate  efficiency  at  6  MV  is 
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more  than  a  factor  of  2  higher  near  180°  than  at  0°  (dose  at 
exactly  180°  decreases  due  to  self-absorption  in  the  rod). 
This  suggests  that,  to  maximize  the  extracted  dose,  future 
rod-pinch  diode  experiments  at  or  above  6  MV  should  be 
performed  in  negative  polarity  with  an  appropriately 
designed  anode  to  minimize  self-absorption.  [8]  Also 
shown  in  Fig.  3  are  power-law  fits  to  the  dose-rate  at 
various  angles  of  the  form  A  I  V71.  These  fits  apply  to 
voltages  up  to  6  MV  and  are  used  to  model  the  Asterix 
experiments  where  peak  voltages  ranged  between  4-  and 
6-MV. 


Figure  3.  Dependence  of  the  dose-rate  efficiency  on 
voltage  for  various  angles.  Open  symbols  are  from 
simulations  and  dashed  lines  are  power-law  fits  for 
voltages  up  to  6  MV. 

Calculated  and  measured  dose-rates  for  a  typical 
Asterix  shot  are  compared  in  Fig.  4.  The  calculated  dose- 
rates  are  obtained  using  the  measured  current  and  voltage 
waveforms  (shown  in  Fig.  4 A)  in  the  power-law  fits  from 
the  simulations.  The  diode  voltage  is  obtained  by 
applying  a  rather  large  inductive  correction  (L=1030  nH) 
to  an  electrolytic  resistive  divider  located  upstream  of  the 
diode. [3]  As  a  result,  there  is  a  rather  large  uncertainty  in 
the  voltage  waveform.  The  current  is  an  average  of  three 
B-dot  loops  located  near  the  diode  load.  The  variation  of 
the  individual  B-dots  from  this  average  is  small  indicating 
a  less  than  ±5%  uncertainty  on  the  current.  The 
magnitudes  and  shapes  of  the  dose-rates  calculated  from 
the  LSP/CYLTRAN  model  are  compared  with  the 
measured  dose-rates  at  0°,  20°,  40°,  60°,  and  80°  in 
Figs.  4B-4F.  A  comparison  of  the  angular  variation  of  the 
total  calculated  dose  (i.e.  the  time-integral  of  the  dose- 
rate)  and  measured  the  measured  TLD  dose  are  shown  in 
Fig.  5.  The  calculated  doses  and  dose-rates  are  in 
agreement  with  the  measurements  to  within  the 
experimental  precision.  Similar  agreement  is  obtained  for 
a  shot  that  used  1  -mm  diameter  tapered  tungsten  rod. 

A  comparison  of  the  calculated  dose-rate  efficiencies 
from  1-mm  (Fig.  6A)  and  2-mm  (Fig.  6B)  diameter 
tapered  and  blunt  rods  are  shown  in  Fig.  6.  For  a  given 
rod  diameter,  this  figure  shows  that  the  0°  dose-rate 


calculated  dose-rates  for  a  2-mm  diameter,  tapered 
tungsten  rod.  A)  The  voltage  and  current  waveforms. 
B)  0°  dose-rates.  C)  20°  dose-rates.  D)  40°  dose-rates. 
E)  60°  dose-rates.  F)  80°  dose-rates. 

efficiency  is  higher  for  tapered  rods  than  blunt  rods.  Near 
180°,  the  dose-rate  efficiency  from  blunt  rods  is 
significantly  higher  than  the  tapered  rods.  The  dose-rate 
efficiency  at  10  MV  is  also  seen  to  be  more  than  ten  times 
higher  near  180°  than  at  0°.  This  suggests  that  blunt  rods 
are  better  able  to  take  advantage  of  the  electron  angles  of 
incidence  as  the  voltage  increases.  The  calculations  at 
10  MV  indicate  a  current  of  Fr=210  kA  for  a  ratio  of  the 
cathode-to-anode  radii  of  rc/rA=ll.  Therefore,  a  dose  of 
-800  Rad  (LiF)  near  180°  is  possible  1  m  from  a  2-mm 
diameter,  blunt  rod  in  a  50-ns  full-width  at  half  maximum 
radiation  pulse.  Further  increases  in  the  dose  are  possible 
if  the  ion  current  fraction  (-40%  at  10  MV)  can  be 


Figure  5.  Comparison  of  calculated  and  measured  angular 
dose  variation  for  a  2-mm  diameter,  tapered  tungsten  rod. 
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Figure  6.  Comparison  of  the  calculated  0°  and  180°  dose- 
rate  efficiencies  for  blunt  and  tapered  rods.  A)  1-mm 
diameter  rods.  B)  2-mm  diameter  rods. 

reduced  by  eliminating  protons.  When  protons  are 
replaced  by  tungsten  ions,  the  ion  current  fraction  is 
reduced  to  less  than  5%  and  the  dose  increases  by  50%. 
One  method  for  eliminating  protons  is  to  heat  the  tungsten 
anode  to  about  2000°  K  prior  to  the  shot  to  drive  off 
impurities.  [9] 

The  simulations  indicate  very  little  difference  in  the 
dose-rate  efficiency  from  1-  and  2-mm  diameter  blunt 
rods  at  180°.  This  can  be  understood  by  comparing  the 
current  flow  patterns  from  1-  and  2-mm  diameter  blunt 
rod  pinch  diodes  at  10  MV  shown  in  Fig.  7.  From  this 
figure  it  is  seen  that  electrons  E*B  drift  down  the  length 
of  the  rod  along  the  virtual  cathode  surface  defined  by 
|E|~0.  The  strong  self-magnetic  field  causes  the  current 
flow  to  bend  around  the  end  of  the  rod  and  the  majority  of 
the  current  flow  enters  the  rod  perpendicular  to  the  end  of 
the  rod  for  both  the  1-mm  diameter  (Fig.  7 A)  and  the  2- 
mm  diameter  (Fig.  7B)  blunt  rods.  In  both  cases  the 
majority  of  the  electron  beam  couples  to  the  rod  axially 
and  not  radially  so  that  the  x-ray  production  efficiency  is 
independent  of  the  rod  diameter.  Side  viewing  x-ray 
images  indicate  that  a  significant  fraction  of  the  electrons 
are  concentrated  at  the  end  of  the  blunt  rod.  [10] 


120  kA  *  & 


Figure  7.  Current  flow  patterns  from  LSP  simulations  of 
a  10  MV,  blunt- anode,  rod-pinch  diode  with  rc/rA=l  1.  A) 
1-mm  diameter  rod.  B)  2-mm  diameter  rod. 
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